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Abstract 

A  three-dimensional,  single-phase,  multi-component  mathematical  model  has  been  developed  for  a  liquid-fed  direct  methanol  fuel  cell  (DMFC). 
The  traditional  continuity,  momentum,  and  species  conservation  equations  are  coupled  with  electrochemical  kinetics  in  both  the  anode  and  cathode 
catalyst  layer.  At  the  anode  side,  the  liquid  phase  is  considered,  and  at  the  cathode  side  only  the  gas  phase  is  considered.  Methanol  crossover  due 
to  both  diffusion  and  electro-osmotic  drag  from  the  anode  to  the  cathode  is  taken  into  consideration  and  the  effect  is  incorporated  into  the  model 
using  a  mixed-potential  at  the  cathode.  A  finite-volume-based  CFD  technique  is  used  to  develop  the  in-house  numerical  code  and  the  code  is 
successfully  used  to  simulate  the  fuel  cell  performance  as  well  as  the  multi-component  behavior  in  a  DMFC.  The  modeling  results  of  polarization 
curves  compare  well  with  our  experimental  data.  Subsequently,  the  model  is  used  to  study  the  effects  of  methanol  crossover,  the  effects  of  porosities 
of  the  diffusion  layer  and  the  catalyst  layer,  the  effects  of  methanol  flow  rates,  and  the  effects  of  the  channel  shoulder  widths. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

The  direct  methanol  fuel  cell  (DMFC)  is  considered  a  highly 
promising  power  source  with  its  important  attributes  of  quick 
refueling,  low  temperature  and  pressure,  low  methanol  cost,  no 
liquid  electrolyte,  and  compact  cell  design,  etc.  Efforts  in  devel¬ 
oping  mathematical  models  for  DMFCs  have  been  limited  until 
recent  years.  Baxter  et  al.  [1]  developed  a  one-dimensional  and 
single-phase  mathematical  model  for  a  liquid-fed  DMFC,  focus¬ 
ing  on  the  anode  catalyst  layer.  Dhole  et  al.  [2]  presented  a  one¬ 
dimensional  model  for  the  vapor-fed  DMFC  and  the  crossover 
effects  were  studied.  Scott  et  al.  [3,4]  developed  several  simpli¬ 
fied  single-phase  models  to  study  transport  and  electrochemical 
processes  in  a  liquid-fed  DMFC  and  showed  that  the  cell  per¬ 
formance  is  limited  by  the  slow  diffusion  of  methanol  in  liquid. 
Scott  and  co-workers  [5]  also  developed  a  model  to  predict  the 
effective  methanol  concentration  at  the  catalyst  surface.  This 
model,  together  with  an  empirical  model  of  the  open  circuit  volt¬ 
age  and  the  cathode  over-potential  model,  was  used  to  predict  the 
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overall  cell  voltage  and  current  density.  Cruickshank  and  Scott 
[6]  presented  a  simple  model  to  describe  methanol  permeation. 
Shukla  et  al.  [7]  developed  a  one-dimensional  model  to  com¬ 
pare  the  performance  of  a  solid-polymer  electrolyte  DMFC  with 
aqueous  methanol  and  methanol  mixed  with  air  at  the  anode  side. 
Kulikovsky  et  al.  [8]  simulated  a  vapor-fed  DMFC  with  a  two- 
dimensional  model  and  compared  the  detailed  current  density 
distributions  in  the  backing,  catalyst  layer,  and  membrane  sepa¬ 
rator  between  a  conventional  and  a  new  current  collector.  Wang 
and  Wang  [9]  developed  two-dimensional,  two-phase  and  multi- 
component  model  for  a  liquid-fed  DMFC.  In  this  model,  they 
considered  the  mixed  potential  effects  of  methanol  oxidation 
due  to  methanol  crossover,  but  the  catalyst  layers  were  treated  as 
infinitely  thin  interfaces.  Birgersson  et  al.  [10]  developed  a  two- 
dimensional  and  single-phase  model  for  the  anode  of  a  DMFC. 
In  their  model,  the  catalyst  layer  was  taken  into  consideration 
as  boundary  conditions  via  parameter  adaptation.  Divisek  et  al. 
[11]  accounted  for  the  influence  of  both  the  methanol  and  oxy¬ 
gen  reaction  kinetics  and  their  dependency  on  the  two-phase 
mass  flow  in  the  catalyst  and  diffusion  layers  by  a  vapor-liquid 
model  for  the  DMFC.  In  this  work,  we  have  developed  a  three- 
dimensional,  single-phase  and  multi-component  model  for  the 
entire  DMFC,  including  both  the  anode  and  cathode  channels, 
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Nomenclature 

a  water  activity 

ai!jef  reference  exchange  current  density  times  specific 

area  at  anode  (Am-3) 

ai£ef  reference  exchange  current  density  time  specific 
area  at  cathode  (Am-3) 

C  molar  concentration  (mol  m-3) 

D  diffusion  coefficient  (m2  s-1) 

E  voltage  (V) 

F  Faraday’s  constant  (C  mol  - 1 ) 

hc  channel  depth  (z  direction  width)  (m) 

I  current  density  (A  m-2) 

j  transfer  current  density  (A  m-3) 

k^  electrokinetic  permeability  of  membrane  (m2) 

kh  hydraulic  permeability  (m2) 

&air  permeability  to  air  of  gas  diffuser  (m2) 
kw  permeability  to  water  of  gas  diffuser  (m2) 

Lc  channel  length  (y  direction  width)  (m) 

M  molecular  weight  (g  mol- 1 ) 

N  mole  flux  (mol  m-2  s- 1 ) 

P  pressure  (N  m-3) 

Q  flow  rate  (m3  s  - 1 ) 

rs  porous  media  correction  factor 

R  resistance  (Q) 

Rc  gas  constant  (J  mol- 1  K- 1 ) 

S  source  term 

fa  diffusion  layer  thickness  (m) 

tc  catalyst  layer  thickness  (m) 

tm  membrane  thickness  (m) 

T  temperature  (K) 

V  velocity  (ms-1) 

wc  half-channel  width  (m) 

ws  collector  plate  width  (m) 

X  mole  fraction 

z  charge  of  fixed  (sulfonate)  sites 

Greek  symbols 
a  transfer  coefficient 

s  porosity 

y  reaction  order 

A  electro-osmotic  drag  coefficient 

li  viscosity  (kgm-1  s-1) 

p  density  (kgm-3) 

o  conductivity  (l/(£2  m)) 

\f/  water  content  in  the  membrane 


Subscripts  and  superscripts 


a 

anode,  air 

c 

cathode 

ca 

catalyst  layer 

cross 

methanol  crossover 

CH 

methanol 

d 

diffusion 

eff 

effective 

g  gas 

H+  proton 

i  ith  layer 

k  klh  component 

1  liquid  water 

m  membrane 

O  oxygen 

ref  reference 

w  water 


both  diffusion  layers,  and  both  catalyst  layers  and  the  membrane. 
The  model  for  the  anode  and  cathode  sides  is  coupled  by  current 
density.  The  modeling  results  of  polarization  curves  compared 
well  with  our  experimental  data  and  this  model  is  used  to  study 
the  effects  of  methanol  crossover  and  various  parameter  effects. 

2.  Mathematical  model 

2.7.  Model  description 

The  modeling  domain  for  the  DMFC  consists  of  the  anode 
side,  the  cathode  side,  and  the  membrane  as  shown  in  Fig.  1. 
Both  the  anode  and  cathode  sides  have  the  same  structure  and 
can  be  sub-divided  into  the  gas  channel,  the  diffusion  layer,  and 
the  catalyst  layer.  In  Fig.  1,  /zc  is  the  gas  channel  height;  fa,  tc , 
tm  the  thickness  of  the  diffusion  layer,  the  catalyst  layer,  and  the 
membrane,  respectively;  wc,ws,  wt  are  the  half- widths  of  the  gas 
channel,  the  collector  plate,  and  the  diffusion  layer,  respectively; 
and  Lc  is  the  length  along  the  channel.  The  reactants  enter  the 
gas  channels  from  the  surface  at  y  =  0.  The  detailed  geometric 
parameters  are  given  in  Table  1 . 

Two  components,  methanol  and  water,  were  considered  at  the 
anode  side,  and  air  was  considered  at  the  cathode  side. 


Fig.  1.  Schematic  of  the  modeling  domain. 
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Table  1 


Geometric  parameters  used  in  the  model 


Channel  length  (Lc) 

Half-channel  width  (wc) 

Channel  height  ( hc ) 

Channel  shoulder  width  (ws) 

Diffusion  layer  thickness  (p) 

Catalyst  layer  thickness  (tc) 

Membrane  thickness  (Nafion®117)  (tm) 

6.5  x  10-2  m 
5.0  x  10-4  m 
8  x  10-4  m 
1.0  x  10-3  m 
1.4  x  10~4  m 
3  x  10-5  m 
1.8  x  10-4  m 

Table  2 

rs  and  Sj  in  momentum  equations 

G 

Si 

Channel  1.0 

0 

Diffuser  2.25  (  ^  —  1  j 

Catalyst  2.25 (£  -  l)2 

-e^V 
-*c£v  + 

tc*+zF ^ 

2.2.  Model  assumptions 

In  the  model  development,  the  following  assumptions  are 
used:  fluids  are  incompressible;  all  the  flows  are  laminar  and  all 
the  processes  are  steady  state  as  well  as  isothermal;  the  effects 
of  carbon  dioxide  generated  in  the  anode  side  are  negligible; 
methanol  crossed  over  from  the  anode  to  the  cathode  is  com¬ 
pletely  oxidized  at  the  interface  between  the  membrane  and  the 
cathode  catalyst  layer;  and  the  membrane  is  impermeable  to 
gases  and  is  fully  hydrated. 


Table  3 


Anode  and  cathode  source  term  (sf)  in  species  equations 


Anode  (source  term  sk) 

Cathode  (source  term  sk) 

Methanol 

Water 

Proton 

Oxygen 

Water 

Proton 

Channel  layers 

0 

0 

0 

0 

0 

0 

Diffusion  layers 

0 

0 

0 

0 

0 

0 

Catalyst  layers 

6 F  la 

Mw  • 

6  F  1& 

~p  ja 

Mq  • 

4  F  1 c 

Mw  • 
IF  J c 

~P  jc 

2.3.  Model  equations 

On  both  the  anode  and  cathode  sides,  the  governing  equa¬ 
tions  for  the  flow  field  and  species  concentration  are  in  the  same 
form.  The  appropriate  model  equations  in  vector  forms  are  listed 
below. 


Continuity  equation 

< 

II 

o 

(1) 

Momentum  equation 

pV  ■  VV  =  -VP  +  r8n  V2V  +  pSj 

(2) 

In  the  momentum  equations  and  Table  2,  s&  and  sc  are  the 
porosities  of  the  diffusion  layer  and  catalyst  layer,  respectively, 
r£  the  porous  media  correction  factor,  k  the  permeability  of  water 
(kw)  in  the  anode  and  the  permeability  of  air  (ka\T)  in  the  cath¬ 
ode,  kef)  the  electrokinetic  permeability  of  the  membrane,  kh  the 
hydraulic  permeability  of  the  membrane,  z  the  charge  number  of 


Table  4 

Physical  parameters  and  basic  operation  conditions 


T=  343  K 

P=  1.013  x  105  N m-2 

ginid  _  iq-5  m3  S-1 

Qlnlet  =  2  x  10-5  m3  s-1 
Aglet  =  0.21 


C^et  =  1000  mol  m-3 

C“f  =  2000  mol  m“3 

Cgf  =  0.472  mol  m“3 

[7] 

D0  =  1.22  x  10“10m2s“1 

Dch  =  2.8  x  10-9  m2  s_1 

[4] 

Cell  temperature 
Pressure 

Anode  reactant  flow  rate 
Cathode  reactant  flow  rate 
Inlet  oxygen  mole  fraction 
Inlet  methanol  concentration 
Methanol  reference  concentration 
Oxygen  reference  concentration 
Diffusion  coefficient  of  oxygen  in  air 
Diffusion  coefficient  of  methanol  in  water 
Porosity  of  diffusion  layer 
Porosity  of  catalyst  layer 
Porosity  of  membrane 
Permeability  to  air  in  the  gas  diffuser 
Permeability  to  water  in  the  gas  diffuser 
Air  density 
Water  density 
Air  viscosity 

Electrokinetic  permeability  of  the  membrane 
Hydraulic  permeability  of  the  membrane 
Gas  constant 

Reference  exchange  current  density  times  specific  area  at  anode 

Reference  exchange  current  density  times  specific  area  at  cathode 

Anode  reaction  order 

Cathode  reaction  order 

Electro-osmotic  drag  coefficient  of  water 

Anode  transfer  coefficient 

Cathode  transfer  coefficient 

Charge  of  fixed  (sulfonate)  sites 


£d  =  0.6 

£c  =0.4 

£m  =  0.28  [16] 

kair  =  1.76  x  10-11  m2 

kw  =  1.0  x  10-11  m2 

Pair  =  1.0  kg  m-3 

pw  =  1000  kg  m-3 

fig  =  2.05  x  10-5  kgm-1  s_1 


4  =  7.18  x  10"2°m2 

[16] 

4  =  1-8  x  10"l8m2 

[16] 

Sc  =  8.314Jmol-1  K-1 

ai',cl  =  1.0  x  106  Am-3 

[4] 

ai^ef  =  200  Am-3 

[7] 

Ya  =  1.0 

[7] 

Yc  =  1.0 

[7] 

A.w  —  2.5 

[1] 

aa  =  0.5 

[7] 

in 

o 

ii 

o 

z—  —  1 
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the  fixed  sites,  cH+  the  concentration  of  protons  that  is  taken  to 
be  the  concentration  of  the  fixed  charge  and  0  is  the  membrane 
phase  potential. 


Species  equations 

p  V  ■  VXk  =  pDf  V2Xk  +  Sk 

(3) 

=  i 

(4) 

k 


In  Eqs.  (3)  and  (4)  as  well  as  Table  3,  k  represents  methanol 
or  water  in  the  anode  side,  and  oxygen,  nitrogen,  or  water  vapor 
in  the  cathode  side,  Xk  the  mole  fraction  of  species  k,  Sk  the  mass 
generation  rate  for  species  k  per  unit  volume,  Dff  the  effective 
diffusion  coefficient  of  the  kth  component  and  ya  and  jc  are  the 
volumetric  current  density  of  anode  and  cathode,  respectively. 
The  effective  diffusivity  is  given  by: 

Df  =  DksL5  (5) 

where  s  is  the  porosity  of  porous  media  and  Dk  is  the  diffusivity 
of  the  kth  component. 


2.3.1.  Electrochemical  kinetics  in  catalyst  layers 

The  Tafel  equation  is  used  to  determine  the  transfer  current 
densities  jc  and  ja: 


•ref 

Ja  =  aia 


XCH 

vref 

ACH 


Fa 

exp 


a> 

CD 

(0 

o 

> 


■ 

1.0M  experiment 

1.0M  model 

o 

2.0M  experiment 

2.0M  model 

A 

0.5M  experiment 

0.5M  model 

A 


0.3 


Curent  density  (A/cm2) 


0.4 


o  "■ 


~i — 

0.5 


Fig.  2.  Comparison  of  modeling  results  with  experimental  data  at  different 
methanol  feed  concentrations.  The  cell  temperature  70  °C,  cathode  humidi¬ 
fication  temperature  70  °C,  methanol  flow  rate  6  ml  min-1,  and  air  flow  rate 
1200  seem. 


where  crjjf  is  the  reference  ionic  conductivity  at  303  K. 

The  resistance  of  the  membrane  is  defined  by: 

ftm  1  fin 

Rm=  —  dz=—  (11) 

JO  ^fin  Cfin 

where  tm  is  the  thickness  of  the  membrane  and  am  is  the  mem¬ 

brane  conductivity. 


Xp 

vref 

Ao 


Yc 

exp 


where  ai!jef  and  ai£ef  are  the  reference  exchange  currents  density 
multiplied  by  the  specific  area  at  the  anode  and  cathode  sides, 
respectively,  ya  and  yc  the  anode  and  cathode  reaction  order,  rja 
and  ijc  the  anode  and  cathode  over-potential  and  aa  and  ac  are 
the  anode  and  cathode  transfer  coefficients,  respectively. 


2.3.3.  Methanol  crossover 

The  methanol  transport  is  treated  similarly  to  the  water  trans¬ 
port.  In  the  membrane,  the  methanol  transfer  is  caused  by  dif¬ 
fusion  and  electro-osmotic  drag,  and  the  methanol  flux  is  given 
by  [4]: 


^CH 


k-CH  I 
F 


—  c1-5  nm 
£m  UC\  I 


dCgH 

d  z 


(12) 


2.3.2.  Cell  potential  and  current  density 
The  cell  voltage  is  calculated  by: 

^cell  =  E0  fia  T  fic  kRm  (8) 


where  Eo  is  the  open  circuit  voltage,  Rm  the  membrane  resistance 
and  I  is  the  average  current  density  of  the  cell  and  is  given  by: 

1=  Tiadz  (9) 

Jo 

where  tc  is  the  catalyst  layer  thickness. 

The  membrane  is  in  direct  contact  with  liquid  water  at  the 
anode  side  and  water  is  produced  at  the  cathode  side.  Besides, 
water  is  also  supplied  to  the  cathode  side  by  electro-osmosis. 
Thus,  it  is  a  reasonably  to  assume  the  membrane  is  always  fully 
hydrated  and  for  a  fully  hydrated  membrane  its  ionic  conductiv¬ 
ity  can  be  determined  by  [12]: 


crm(T)  =  aTf  ex  p 


(10) 


0.0000  0.0002  0.0004  0.0006  0.0008  0.0010 

z-direction  distance  from  channel  (m) 


Fig.  3.  Average  methanol  mole  fraction  profiles  at  different  current  densities. 
The  cell  temperature  70  °C,  methanol  feed  concentration  0.5  M,  methanol  flow 
rate  6  ml  min-1 ,  and  air  flow  rate  1200  seem. 
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where  C^H  is  the  methanol  concentration  in  the  membrane.  The 
electro-osmotic  drag  coefficient  for  methanol,  Ach,  is  defined 
as  the  number  of  methanol  molecules  dragged  by  each  proton 
through  the  membrane,  (CH30H/H+),  and  it  is  given  by  [13]: 

^CH  =  ^CHlac/m^w  (13) 

where  Xqr  |ac/m  is  the  methanol  mole  fraction  at  the  interface 
between  the  anode  catalyst  layer  and  the  membrane. 

The  diffusion  coefficient  of  methanol  in  the  membrane  (T)^H) 
is  given  by  [4] : 

D£h  =  4.012  x  10“13  exp(0.0243 12T)  (m2  s_1)  (14) 

2.4.  Solution  method 

The  velocity  fields  at  the  anode  and  cathode  sides  are  deter¬ 
mined  first.  For  a  given  value  of  the  anode  over-potential  ^a, 
the  species  concentration  and  the  average  current  density  at  the 


anode  side  are  determined.  Then  the  flux  of  methanol  crossover 
can  be  calculated  and  the  pseudo  current  density  due  to  methanol 
crossover  is  determined  by: 

/cross  =  —ffl  (15) 

6  F 

Using  species  equations  and  the  Tafel  equation  in  the  cathode 
side  iteratively  [9] : 

('  + exp ("lU  m 

the  cathode  over-potential  can  be  determined  and  the  cell  voltage 
is  calculated  using  Eq.  (8).  Then  another  value  of  r/a  is  given 
and  the  process  is  repeated  until  the  entire  polarization  curve  is 
obtained. 


(b)  Cell  current  density  (A/cm2) 


(c)  Cell  current  density  (A/cm2) 

Fig.  4.  (a)  Methanol  concentration  vs.  cell  current  density  at  the  anode/membrane  interface.  The  cell  temperature  is  70  °C,  methanol  flow  rate  6  ml  min-1  and  air 
flow  rate  1200  seem,  (b)  Equivalent  current  density  vs.  cell  current  density.  The  cell  temperature  70  °C,  methanol  flow  rate  6  ml  min-1  and  air  flow  rate  1200  seem, 
(c)  Equivalent  current  density  of  methanol  crossover  vs.  cell  current  density  at  different  methanol  concentrations.  The  cell  temperature  70  °C,  methanol  flow  rate 
6  ml  min-1 ,  and  air  flow  rate  1200  seem. 
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3.  Modeling  results  and  discussions 

3.1.  Comparison  with  experimental  data 

The  geometric  parameters  used  in  the  model  are  the  same 
as  those  used  in  our  previous  experimental  studies  [14]  and  are 
listed  in  Table  1 .  The  physical  parameters  and  basic  operating 
conditions  are  tabulated  in  Table  4.  Typical  modeling  results  of 
cell  polarization  curves  at  different  methanol  concentrations  are 
shown  in  Fig.  2  and  compared  with  the  experimental  results  [  14] . 
The  cell  temperature  is  70  °C,  cathode  humidification  tempera¬ 
ture  is  70  °C,  methanol  flow  rate  is  6  ml  min-1 ,  and  air  flow  rate 
is  1200  seem.  It  is  seen  from  Fig.  2  that  the  three-dimensional 
modeling  results  show  good  agreement  with  the  experimental 
data. 

Fig.  3  shows  average  methanol  fraction  profiles  at  three  dif¬ 
ferent  cell  current  densities.  The  modeling  result  suggests  that 
the  gas  diffusion  layer  contributes  a  substantial  fraction  of  the 
total  mass  transfer  resistance,  which  indicates  the  need  of  work 
and  improvement  in  the  gas  diffusion  layer.  This  result  is  in 
agreement  with  those  obtained  by  Scott  et  al.  [3,4]  that  the  cell 
performance  is  limited  by  the  slow  diffusion  of  methanol  in  liq¬ 
uid  water. 

3.2.  Effect  of  methanol  crossover 

Fig.  4(a)  shows  the  methanol  concentration  contribution  at 
different  current  densities  on  the  interface  between  the  anode 
catalyst  layer  and  the  membrane.  Fig.  4(b)  shows  the  cell  current 
density  versus  the  equivalent  crossover  current  densities,  which 
includes  both  the  electro-osmotic  drag  flux  and  the  diffusion 
flux.  Fig.  4(c)  shows  the  modeling  results  of  equivalent  crossover 
current  densities  at  different  cell  current  densities. 

Generally  speaking,  the  amount  of  methanol  crossover 
decreases  with  current  density,  and  this  trend  is  in  good  agree¬ 
ment  with  those  reported  in  the  literatures  [2,15].  As  mentioned 
earlier,  methanol  crossover  is  mainly  due  to  the  diffusion  and 
the  electro-osmosis.  The  diffusion  flux  is  proportional  to  the 
methanol  concentration  gradient  across  the  membrane.  At  the 
interface  between  the  anode  catalyst  layer  and  the  membrane 
methanol  concentration  decreases  with  current  density.  Thus  it 
is  easy  to  understand  that  the  diffusive  methanol  flux  decreases 
with  current  density.  As  shown  in  Eqs.  (12)  and  (13),  the 
methanol  crossover  flux  due  to  electro-osmosis  is  proportional 
to  both  the  current  density  and  the  methanol  concentration  at 
this  interface  since  the  methanol  concentration  is  assumed  to  be 
zero  at  the  membrane  interface  with  the  cathode  catalyst  layer. 
Therefore,  lower  concentration  of  methanol  at  this  interface  may 
also  leads  to  lower  electro-osmosis  flux.  As  the  current  den¬ 
sity  increases  to  such  a  value  that  the  methanol  is  depleted  in 
the  anode  catalyst  layer,  the  methanol  crossover  rate  becomes 
zero.  The  zero  diffusion  flux  is  self-evident  since  the  concen¬ 
tration  gradient  is  zero.  When  the  methanol  concentration  at  the 
interface  between  the  anode  catalyst  layer  and  the  membrane 
becomes  zero,  no  matter  how  much  water  is  transferred  due  to 
electro-osmosis,  no  methanol  is  present  in  water,  thus  methanol 
transfer  due  to  electro-osmosis  is  also  zero. 


Fig.  5.  Modeling  results  of  the  effect  of  methanol  crossover  with  air  flow  on 
cell  performance.  The  cell  temperature  is  70  °C,  methanol  concentration  1  M, 
methanol  flow  rate  6  ml  min- 1 ,  and  oxygen  flow  rate  1200  seem. 

In  Fig.  5,  the  polarization  curve  for  a  fuel  cell  without  the 
effects  of  methanol  crossover  is  compared  with  that  with  the 
effect  of  methanol  crossover.  It  is  clearly  seen  that,  even  at  a 
very  low  methanol  concentration,  e.g.  1  M,  the  effect  of  methanol 
crossover  is  still  significant. 

3.3.  Effect  of  porosities  of  diffusion  layer  and  catalyst  layer 

Figs.  6-9  show  the  effects  of  porosities  of  the  diffusion  layer 
and  the  catalyst  layer  on  the  cell  performance  and  methanol 
crossover.  The  cell  temperature  is  70  °C,  methanol  concentration 
is  0.5  M,  methanol  flow  rate  is  6  ml  min-1  and  air  flow  rate  is 
1200  seem.  From  these  results,  it  can  be  seen  that  the  porosities 
of  both  the  diffusion  layer  and  the  catalyst  layer  have  significant 


Current  density  (A/cm2) 

Fig.  6.  Modeling  results  of  the  effect  of  porosity  of  catalyst  layer  on  perfor¬ 
mance.  The  cell  temperature  is  70  °C,  methanol  concentration  0.5  M,  methanol 
flow  rate  6  ml  min-1 ,  and  air  flow  rate  1200  seem. 
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Fig.  7.  Modeling  results  of  the  effect  of  porosity  of  catalyst  layers  on 
methanol  crossover.  The  cell  temperature  is  70  °C,  methanol  concentration 
0.5  M,  methanol  flow  rate  6  ml  min- 1 ,  and  air  flow  rate  1200  seem. 


Fig.  9.  Modeling  results  of  the  effect  of  porosity  of  diffusion  layers  on 
methanol  crossover.  The  cell  temperature  is  70  °C,  methanol  concentration 
0.5  M,  methanol  flow  rate  6  ml  min-1 ,  and  air  flow  rate  1200  seem. 


effects  on  cell  performances.  The  performance  and  methanol 
crossover  increase  with  the  porosity.  These  results  again  indicate 
that  the  slow  diffusion  of  methanol  in  liquid  water  is  a  liming 
factor  in  DMFC  performances. 

Fig.  6  shows  the  effects  of  porosity  of  the  anode  catalyst 
layer  on  the  cell  performances.  The  modeling  suggests  that  the 
catalyst  layer  porosity  has  a  significant  influence  on  the  cell 
performance.  Fig.  7  shows  the  effects  of  porosity  of  the  anode 
catalyst  layer  on  methanol  crossover.  Fig.  8  shows  the  effects  of 
porosity  of  the  anode  diffusion  layers  on  the  cell  performances. 
Fig.  9  shows  the  effects  of  porosity  of  the  anode  diffusion  layer 
on  methanol  crossover.  Generally  speaking,  methanol  crossover 
decreases  with  decrease  in  porosity  in  the  anode  diffusion  layer; 
while  the  cell  performance  increases  with  the  increase  in  porosity 
in  both  the  anode  diffusion  layer  and  catalyst  layer. 


3.4.  Effect  of  the  methanol  flow  rate 

Fig.  10  shows  the  polarization  curves  at  different  methanol 
flow  rates.  The  cell  temperature  is  70  °C,  methanol  concentra¬ 
tion  is  0.5  M,  and  air  flow  rate  is  1200  seem.  The  current  density 
increases  with  the  anode  methanol  flow  rate  at  the  methanol 
concentration  of  0.5  M.  When  the  anode  flow  rate  is  higher  than 
7  ml  min- 1 ,  the  cell  current  density  does  not  change  significantly 
with  the  anode  flow  rate.  At  low  flow  rate,  due  to  mass  transfer 
resistance,  the  methanol  concentration  is  too  low  in  the  catalyst 
layer  and  causes  a  lower  current  density,  especially  at  the  down 
stream.  When  the  flow  rate  is  high  enough,  any  further  increase 
in  the  flow  rate  has  no  significant  effect  on  the  methanol  concen¬ 
tration  in  the  channel  and  in  the  catalyst  layer,  thus  no  significant 
effect  on  cell  current  density. 


Fig.  8.  Modeling  results  of  the  effect  of  porosity  of  diffusion  layers  on 
methanol  crossover.  The  cell  temperature  is  70  °C,  methanol  concentration 
0.5  M,  methanol  flow  rate  6  ml  min- 1 ,  and  air  flow  rate  1200  seem. 


Current  density  (A/cm2) 

Fig.  10.  Modeling  results  of  polarization  curves  of  different  methanol  feed  flux. 
The  cell  temperature  is  70  °C,  methanol  concentration  0.5  M,  and  air  flow  rate 
1200  seem. 
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Fig.  11.  Modeling  results  of  the  effect  of  channel  shoulder  width  on  perfor¬ 
mance.  The  cell  temperature  is  70  °C,  methanol  concentration  0.5  M,  methanol 
flow  rate  6  ml  min-1 ,  and  air  flow  rate  1200  seem. 


Fig.  12.  Modeling  results  of  the  effect  of  channel  shoulder  width  on  methanol 
crossover.  The  cell  temperature  is  70  °C,  methanol  concentration  0.5  M, 
methanol  flow  rate  6  ml  min- 1 ,  and  air  flow  rate  1200  seem. 

3.5.  Effect  of  the  channel  shoulder  width 

Fig.  11  shows  the  effect  of  the  channel  shoulder  width  on 
the  cell  performance.  It  can  be  seen  that  the  cell  performance 
decreases  with  the  channel  shoulder  width,  but  the  effects  are 
not  significant.  Fig.  12  shows  the  effect  of  the  channel  shoulder 
width  on  methanol  crossover.  The  crossover  decreases  with  the 
increase  in  shoulder  width  due  to  the  blockage  effect  of  the 
shoulders. 

4.  Concluding  remarks 

A  three-dimensional,  single-phase,  multi-component  model 
was  developed  for  a  liquid-fed  direct  methanol  fuel  cell  (DMFC). 


The  electrochemical  kinetics  in  both  the  cathode  and  anode  were 
coupled  with  the  traditional  continuity,  momentum,  and  species 
conservation  equations,  and  the  effect  of  methanol  crossover 
was  also  incorporated  in  the  model  using  a  mixed-potential  at 
the  cathode.  A  finite- volume-based  CFD  technique  was  used 
to  develop  the  in-house  numerical  code  and  the  code  was  suc¬ 
cessfully  used  to  simulate  the  fuel  cell  performance  and  the 
multi-component  behavior  in  a  DMFC.  The  model  was  able  to 
predict  polarization  curves  at  various  different  operating  condi¬ 
tions  and  it  was  used  to  study  the  effects  of  methanol  crossover, 
the  effects  of  porosities  of  diffusion  layers  and  catalyst  lay¬ 
ers,  the  effects  of  methanol  flow  rates,  and  the  effects  of  the 
channel  shoulder  widths.  The  simulation  results  showed  that 
the  effect  of  methanol  crossover  is  significant,  even  at  a  very 
low  methanol  concentrations;  the  methanol  crossover  flux,  thus 
the  effect  on  cell  performance,  decrease  with  increasing  current 
density;  both  the  porosities  of  the  anode  catalyst  layer  and  the 
diffusion  layer  have  significant  effect  on  methanol  crossover  and 
the  cell  performance;  the  effect  of  flow  field  shoulder  width  is  not 
significant. 

It  is  known  that  both  gas  and  liquid  phases  exist  in  the  anode 
and  cathode  in  a  real  DMFC  and  an  accurate  model  needs  to 
take  such  two-phase  phenomena  into  consideration.  Develop¬ 
ment  of  such  a  multi-dimensional  and  multi-phase  flow  model 
is  underway. 
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